APS/123-QED 



Current-induced dynamics in non-collinear dual spin-valves 

P. Balaz, 1 ^ M. Gmitra, 2 and J. Barnas 1 ' 3 

department of Physics, Adam Mickiewicz University, Umultowska 85, 61-614 Poznan, Poland 
2 Institute of Phys., P. J. Safdrik University, Park Angehnum 9, 040 01 Kosice, Slovak Republic 
3 Institute of Molecular Physics, Polish Academy of Sciences Smoluchowskiego 17, 60-179 Poznan, Poland 

(Dated: September 25, 2009) 

Spin-transfer torque and current induced spin dynamics in spin-valve nanopillars with the free 
magnetic layer located between two magnetic films of fixed magnetic moments is considered theo- 
retically. The spin-transfer torque in the limit of diffusive spin transport is calculated as a function 
of magnetic configuration. It is shown that non-collinear magnetic configuration of the outermost 
magnetic layers has a strong influence on the spin torque and spin dynamics of the central free 
layer. Employing macrospin simulations we make some predictions on the free layer spin dynamics 
in spin valves composed of various magnetic layers. We also present a formula for critical current 
in non-collinear magnetic configurations, which shows that the magnitude of critical current can be 
several times smaller than that in typical single spin valves. 

PACS numbers: 



I. INTRODUCTION 

Magnetic multilayers are ones of the most relevant ele- 
ments in the development of cutting-edge technologies. If 
current flows along the axis of a metallic hybrid nanopil- 
lar structure, a spin accumulation builds up in the vicin- 
ity of the normal- mctal/ferromagnet interface. More- 
over, the current produces spin-transfer torque (STT) 
which acts on magnetic moments of the ferromagnetic 
layers^. As a consequence, magnetization of a particu- 
lar layer may be driven to an oscillation mode 3 -, or can 
be switched between possible stable states^. However, 
current-induced control of magnetic moments requires a 
rather high current density. Therefore, reduction of the 
critical current density remains the most challenging re- 
quirement from the point of view of possible applications. 
Since the critical currents are related to STT, which in 
the diffusive transport regime is proportional to spin ac- 
cumulation at the normal-metal/ferromagnet interface, 
one may alternatively rise a question of possible ways to 
enhance the spin accumulation. 

One of the possibilities to decrease the critical current 
density in metallic structures is the geometry proposed 
by L. Berger— , in which the free magnetic layer is located 
between two pinned magnetic layers of opposite magne- 
tizations. Indeed, in such a dual spin valve (DSV) geom- 
etry both interfaces of the free layer can generate STT, 
and this may decrease the critical current several timcs^. 

To examine the influence of an additional magnetic 
layer on the free layer's spin dynamics one needs to find 
first the STT acting on the central layer for arbitrary di- 
rection of its magnetization vector. To do this, we employ 
the spin-dependent diffusive transport approach^, based 
on Valet-Fert description^. In this paper we present a 
comprehensive survey of STT-induced effects in DSVs 
for generally non-collinear magnetic configuration of the 
outermost magnetic films. We also study asymmetric 
exchange-biased DSV, in which magnetic moment of one 
of the outer magnetic layers is fixed to an antifcrromag- 



netic layer due to exchange anisotropy. We show that a 
non-standard (wavy-like) angular dependence of STT - 
originally predicted for asymmetric spin valves^ - can 
also occur in DSV geometry. Such a non-standard an- 
gular dependence of the torque is of some importance 
from the application point of view, as it allows to induce 
steady-state precessional modes without external mag- 
netic field&i^. Furthermore, we examine current-induced 
spin dynamics within the macrospin model, and derive a 
formula for critical currents which destabilize trivial fixed 
points of the central spin's dynamics in non-collinear con- 
figurations of the outermost magnetizations. 

In section [TT1 we describe the model and introduce ba- 
sic formula for STT and spin dynamics in DSVs. Then, 
in section IIIII we study dynamics in a symmetric DSV 
geometry, where we compare STT and spin dynamics in 
DSVs and single spin valves (SSVs). Finally, in section 
IIVI we analyze an exchange-biased DSV structure with 
noncollincar magnetic moments of the outermost mag- 
netic layers. Finally, we conclude in section fVl 



II. MODEL 

We consider a multilayer nanopillar structure, 
Fl/Nl/Fc/Nr /Fr , consisting of three ferromagnetic (F) 
layers separated by normal-metal (N) layers; see Fig. Q] 
Spin moment of the central layer, Fc, is free to rotate, 
while the right Fr and left Fl ferromagnetic layers are 
much thicker and their net spin moments are assumed to 
be fixed for current densities of interest. Fixing of these 
moments can be achieved either by sufficiently strong 
coercieve fields, or by exchange anisotropy at interfaces 
with antiferromagnetic materials. 

In the Landau-Lifshitz-Gilbert phenomenological de- 
scription, the dynamics of a unit vector along the net 
spin moment s of the central (free) magnetic layer is de- 
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FIG. 1: Schematic of a dual spin valve with fixed magnetic 
moments of the outer magnetic layers, Fl and Fr, and free 
central magnetic layer, Fc, separated by normal- metal layers 
Nl and Nr. Here, Sl, Sb,, and s are unit vectors along the 
spin moments of the Fl, Fr, and Fc layers, respectively. 



from conduction electrons within a few interfacial atomic 
layers of the ferromagnetic. Thus, the STT acting on the 
central layer Fc can be calculated asr = (h/2) j^), 
where and j ^ are the spin current components per- 
pendicular to magnetic moment of the free layer and 
calculated at the corresponding left and right normal- 
metal/ferromagnet interfaces. The torque consists of in- 
plane Ty (in the plane formed by magnetic moments of 
the two interacting films) and out-of-plane tj_ (normal to 
this plane) parts, r = Ty + tj_, which have the following 
forms 



scribed by the equation 



ds ds _ 

— hasx — = r. 

dt dt 



(1) 



where a is the Gilbert damping parameter. The right- 
hand side represents the torques due to effective magnetic 
field and spin-transfer, 



|7 g |/i s x H cS 



I7gl 
M s d 



(2) 



where 7 g is the gyromagnetic ratio, is the vacuum per- 
meability, M s stands for the saturation magnetization of 
the free layer, and d denotes thickness of the free layer. 
Considering the thin ferromagnetic free layer of ellipti- 
cal cross-section, the effective magnetic field if e g can be 
written as 



H, 



(3) 



and includes applied external magnetic field (H cxt ), uni- 
axial anisotropy field (-ffani)j an d the demagnetization 
field (i?dem)j where e z is the unit vector along the z-axis 
(easy axis) , compare scheme in Fig. [TJ The demagnetiza- 
tion field can be written in the form ffdem = (s-N)M s = 
(HdxSx, HdySy, HdzS z ), where TV is a diagonal demagne- 
tization tensor. 



Spin-transfer torque 



s x ( a L S L - a R Sj. 



Is X 



Isx ( 6 L S L - 6r.5r 



(4a) 
(4b) 



where I is the current density, and Sl and Sr are the 
unit vectors pointing along the fixed net-spins of the Fl 
and Fr layers, respectively. The parameters ol, or, &l, 
and 6r depend, generally, on the magnetic configuration 
and material composition of the system, and have been 
calculated in the diffusive transport limit^. According to 
first-principles calculations of the mixing conductance^ 2 -, 
the out-of-plane torque in metallic structures is about 
two orders of magnitude smaller than the in-plane torque. 
Although this component has a minor influence on crit- 
ical currents, it may influence dynamical regimes, so we 
include it for completeness in our numerical calculations. 

Let us consider now how the STT affects spin 
dynamics of the free magnetic layer. We rewrite 
Eq.([T]) in spherical coordinates (<j),6), which obey s = 
(cos <j> sin 8, sin </> sin 8, cos 8) ; see Fig. Q] Defining unit 
base vectors, = (e 2 xs)/ sin 8 and eg = e^x s, Eq.|T|) 
can be rewritten as 



d_ (<t>\ 1 / sin 
dt V* 9 / 1 + cr 2 V a 



—a sin 
1 



(5) 



Generally, STT acting on a magnetic layer is deter- 
mined by the electron spin angular momentum absorbed 



where the overall torques vg and v^, changing the angles 
8 and 0, respectively, read 



vg = r • eg = - |7 g |/i (H dx - H dy ) cos (j) sin cj) sin 8 + Tg , (6a) 



= T • = -|7g|/x [H cxt + (Hani + H dx cos 2 4> + H dy sin 2 - H dz ) cos 9] sin 8 + t . (6b) 

I 



The first terms in Eqs. ^ describe the torques due to demagnetization and anisotropy fields of the free layer, 
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while Te = Tq+Tq- and = t^+t^ arc the corresponding 
components of the current-induced torque, which origi- 
nate from tii and Tj_. 

As we have already mentioned above, the main con- 
tribution to STT comes from Ty. Assuming now that 
magnetic moment of the left magnetic layer is fixed 
along the z-axis, Sl = (0,0,1), and magnetic moment 
of the right magnetic layer is rotated by an angle Q from 
the z-axis and fixed in the layer's plane (see Fig. [I]), 
Sr = (0, sinfl, cosfi), one finds 

T e — ( a L — a R cos fl)I sin 8 + ap.7 sin Q sin <j) cos 8 , (7a) 
T l — 0,-rI cos 4> sin fl . (7b) 

The component rj consists of two terms. The first one is 
analogous to the term which describes STT in a SSV. 
However, its amplitude is now modulated due to the 
presence of Fr. In turn, the second term in Eq. (|7aj) 
is nonzero only in noncollincar magnetic configurations. 
From Eq. (|7b[) follows that rjj is also nonzero in non- 
collincar configurations, and only if the magnetization 
points out of the layer's plane (</> ^ tt/2). When mag- 
netic moments of the outer magnetic layers are parallel 
(f2 = 0), Tg = (ah — an) I sin 8. For a symmetric DSV, 
dh(0) = or(0), and hence STT acting on s vanishes. On 
the other hand, in the antiparallel configuration of Sl 
and Sj{ (0 = 7r), the maximal spin torque enhancement 

can be achieved, t! = (cil + or)/ sin 8. 

Similar analysis of tj_ leads to the following formulas 
for Tg and t^: 

T e~ = orI cos (j) sin fi , (8a) 
= — (6l — or cos f2)I sin 8 + 6r cos sin ^> sin f2 . (8b) 

Thus, if the outer magnetic moments are collinear, t^ = 
0, while t^J~ reduces to t^ = — (&l — &R)/sin0 for fi = 0, 
and r^" = — (feL + &R.)^sin# for £7 = 7r. Hence, in symmet- 
ric spin valves, where &l = ^r, vanishes in the parallel 
configuration of the outermost magnetic moments and is 
enhanced in the antiparallel configuration. 

In the following sections we investigate STT and its 
effects on critical current and spin dynamics. We start 
from symmetric spin valves in antiparallel magnetic con- 
figuration (£2 = 7r). Then, we proceed with asymmetric 
exchange-biased dual spin valves. 



III. SYMMETRIC DSVS 

Let us consider first symmetric DSVs with antiparallel 
orientation of magnetic moments of the outermost ferro- 
magnetic films: Sl = (0, 0, 1) and Sr = (0, 0,-1). As in- 
dicated by Eqs. (J7J) and (J5J), such a configuration may lead 
to enhancement of STT in comparison to that in SSVs. 
Thus, let us analyze first STT in two types of structures: 
the double spin valve F(20)/Cu(10)/F(8)/Cu(10)/F(20) 



in the antiparallel configuration, and the corresponding 
single spin valve F(20)/Cu(10)/F(8). The numbers in 
brackets correspond to layer thicknesses in nanometers. 

In Fig. [2] we show the angular dependence of STT for 
DSVs and SSVs, when the vector s changes its orien- 
tation, described by angle 8, in the layer plane (<p = 
tt/2). The magnetic layers made of Permalloy, NigoFe2o 
(Fig. [5^), and of Cobalt (Fig. (2h) are considered. Due 
to the additional fixed layer (F R ), STT in DSVs is 
about twice as large as in SSVs, which is consistent 
with Berger's predictions^. Additionally, the angular de- 
pendence of STT acting on the free layer in Co/Cu/Co 
spin valves is more asymmetric than in Py/Cu/Py. This 
asymmetry, however, disappears in Co/Cu/Co/Cu/Co 
DSVs due to superposition of the contributions from both 
fixed magnetic layers to the STT. 
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FIG. 2: Spin transfer torque Tg in symmetric DSVs, 
F(20)/Cu(10)/F(8)/Cu(10)/F(20), in the antiparallel con- 
figuration, Q — 7r, (solid lines), and STT in SSVs, 
F(20)/Cu(10)/F(8) (dashed lines), where (a) F = Permal- 
loy, (b) F = Cobalt. STT is shown in the units of R|I|/|e|, 
and calculated for <j> = 7r/2. The material parameters as in 
Ref[Tl. 

The enhancement of STT in dual spin valves may lead 
to two important improvements of the spin dynamics: 
reduction of the critical current needed to trigger the spin 
dynamics, and decrease of the switching time. The latter 
is defined as the time needed to switch the magnetization 
from one stable position to the opposite one. The fixed 
points of the dynamics of s are given by the equations 
vg = and = 0. If £1 = ir, they are satisfied for 8 = 0, 
and 8 = tt. Employing the 'zero-trace' stability condition 
of the linearized Landau-Lifshitz-Gilbert equation^, we 
find the critical current destabilizing the initial (8 = 0) 
state in the form 



1° 

c, DSV 



where 



afj,oM s d 



dx 



-Ha 



and 



are calculated for 



(9) 

0. Addition- 



ally, we have omitted here the terms resulting from t_l 
because of their small contribution to the critical cur- 
rent. Equation ((9]) is analogous to the expression for 
critical current in SSV—. According to our calculations, 
the critical current in DSVs with F = Cobalt is 6-times 
smaller than in SSVs, as reported by Berger— . However, 
if F = Permalloy, the introduction of a second fixed mag- 
netic layer reduces the critical current only by a factor 
of 2. This difference arises from the dependence of spin 
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accumulation on spin-flip length which is about 10-times 
longer in Cobalt than in Permalloy^. 

To estimate the switching time in DSVs as well as in 
SSVs we employ the single-domain macrospin approx- 
imation to the central layer. Equations fl}, including 
Eq. ((3]) and Eq. (j4]), completely describe the dynamics 
of central layer's spin, s. In our simulations we assumed 
a constant current of density /. The positive current 
(I > 0) is defined for electrons flowing from Fr towards 
Fl (current then flows from Fl towards Fr); opposite 
current is negative. Apart from this, the demagnetiza- 
tion field of the free layer of elliptical cross-section with 
the axes' lengths 130 nm and 60 nm has been assumed, 
while external magnetic field was excluded, -ff C xt = 0. 
For each value of the current density, from Eq. |T]) we 
found evolution of s starting from an initial state until s 
switched to the opposite state (provided such a switch- 
ing was admitted). In all simulation, the spin was ini- 
tially slightly tilted in the layer plane from the orientation 
s = (0,0, 1), assuming 9 = 1° and </>o = n/2. A success- 
ful switching, with the switching time t s , is counted when 
sl(i s ) < —0.99, where sj(f) is the exponentially weighted 
moving average^, sj(i) = t] s z (t) + (1 — r?)sj(i — At), 
At is the integration step, and the weighting parame- 
ter r} = 0.1. The moving average s7 is calculated only 
in time intervals, where s z (t) remains continuously be- 
low the value of —0.9; otherwise s^(t) = s z (t). Fig. [3J 
compares the switching times in DSVs and in the corre- 
sponding SSVs. In both cases shown in Fig. [31 a con- 
siderable reduction of the switching time is observed in 
DSVs. Similarly as for the critical current, the reduction 
of current required for switching in DSVs with Cobalt 
layers is larger from that in DSVs with Permalloy layers. 
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FIG. 3: Switching time in DSVs 

F(20)/Cu(10)/F(8)/Cu(10)/F(20) in the antiparallel mag- 
netic configuration, fl = n, (solid lines), and in SSVs 
F(20)/Cu(10)/F(8) (dashed lines), where (a) F = Permalloy, 
(b) F = Cobalt. The switching time is shown as a function 
of the normalized current density I/Io, with 7o = 10 s Acm -2 . 
The other parameters as in Fig. 2. 



IV. EXCHANGE-BIASED DSV 

Let us consider now an asymmetric exchange-biased 
DSV structure with an antifcrromagnetic layer IrMn ad- 
jacent to the Fr layer in order to pin its magnetic 



moment in a required orientation, i.e. the structure 
Co(20)/Cu(10)/Py(4)/Cu(4)/Co(10)/IrMn(8). The left 
magnetic layer, Fl = Co(20), is assumed to be thick 
enough so its magnetic moment is fixed, Sl = e z . In 
turn, magnetic moment of the right ferromagnetic layer, 
Fr = Co(10), is fixed in the layer plane at a certain angle 
fl with respect to Sl due to the exchange-bias coupling 
to IrMn. 

From Eqs. (O, ®, an d © follows, that in a general 
case (fl ^ 0), the points = and 9 = tt are no more 
solutions of the conditions for fixed points, vg = and 
i>0 = 0, because of the additional terms in STT, which ap- 
pear in non-collincar situations (discussed in section II). 
These additional terms lead to a nontrivial (9-dcpcndence 
of STT, and to a shift of the fixed points out of the 
collincar positions. 

To analyze this effect in more details, let us consider 
first STT assuming s in the layer plane (</> = it/2). Ac- 
cording to Eq. (ffb]) and Eq. ([8b]). and to the fact that the 
parameters b are much smaller than a, the component 
is very small. In Fig. Ufa) we show the second component 
of the torque, namely Tg , as a function of the angle 9 and 
for different values of the angle fl. The configurations 
where Tg — are presented by the contour in the base 
plane of Fig. Ufa). 

In the whole range of the angle fl, two 'trivial' zero 
points are present. Additionally, for small angles close 
to fl = 0, two additional zero points occur. The ap- 
pearance of these additional zero points closely resem- 
bles non-standard wavy-like STT angular dependence, 
which has been already reported in single spin valves 
with magnetic layers made of different material o 6 ' 8 ' 16 . 
Such a ^-dependence dictates that both zero-current 
fixed points (9 = 0,n) become simultaneously stabi- 
lized (destabilized) for positive (negative) current. This 
behavior is of special importance for stabilization of 
the collinear configurations, and for microwave genera- 
tion driven by current only (without external magnetic 
field)£&il. We note, that in contrast to SSVs, the wavy- 
like ^-dependence in exchange-biased DSVs is related 
rather to asymmetric geometry of the multilayer than 
to bulk and interface spin asymmetries. This trend is 
depicted in Fig. 01(b) , where we show variation of STT 
for different thicknesses of Fr at fl = 0. The wavy-like 
torque angular dependence appears for the thickness of 
Fr markedly different from that of Fl . 

Making use of the above described STT calculations, 
extended to arbitrary orientation of s, we performed nu- 
merical simulations of spin dynamics induced by a con- 
stant current in zero external magnetic field (if e xt = 0). 
The sample cross-section was assumed in the form of an 
ellipse with the axes' lengths 130 nm and 60 nm. As be- 
fore, the numerical analysis has been performed within 
the macrospin framework, integrating equation of mo- 
tion, Eq.(fT]). In the simulation we analyzed the long-term 
current-induced spin dynamics started from the initial 
state corresponding to 6>o = 1°, </>o = ir/2. 

As one might expect, the current-induced spin dynam- 
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FIG. 4: (a) STT acting on the central magnetic layer 
in Co(20)/Cu(10)/Py(4)/Cu(4)/Co(10)/IrMn(8) exchange- 
biased DSV as a function of the angle 9, calculated for 
Q = kiv/A, k = 0,1,2,3,4. The contour plot in the base 
plane corresponds to zeros of tq. (b) Wavy- like STT angular 
dependence in exchange biased DSV for = 0, calculated for 
different thicknesses of the Fr layer. STT is shown in the 
units of ft|7|/|e|. The other parameters as in Fig. 2. 



ics of the free layer depends on the angle O, current den- 
sity I, and on current direction. To designate different 
regimes of STT-induced spin dynamics, we constructed 
a dynamical phase diagram as a function of current and 
the angle f2. The diagram shows the average value (s z ) 
of the z component of the free layer net-spin in a sta- 
ble dynamical regime (see Fig. [5^), as well as its disper- 



sion D(s z ) = y (s 2 z ) — (s z ) (see Fig. \5]p). The average 
value (s z ) provides an information on the spin orienta- 
tion, whereas the dispersion distinguishes between static 
states (for which D = 0) and steady precessional regimes 
(where D > 0), in which the z component is involved. 
For each point in the diagram, a separate run from the 
initial biased state </>o = ir/2 and 6q = 1° was performed. 
In the (s z ) diagram, Fig.^a), one can distinguish three 
specific regions. Region (i) covers parameters for which 
a weak dynamics is induced only: s finishes in the equi- 
librium stable point which is very close to s = (0,0,1). 
As the angle f2 increases, STT becomes strong enough to 
cause switching, as observed in the region (ii). The higher 
the angle f2, the smaller is the critical current needed for 
destabilization of the initial state. For smaller f2, the 
current-induced dynamics occurs for currents flowing in 
the opposite direction, see the region (iii). This behav- 
ior is caused by different sign of STT in the initial state. 
From (s z ) we conclude, that neither of the previously 
mentioned stable states is reached. To elucidate the dy- 
namics in region (iii), we constructed the corresponding 
map of D(s z ), Fig.[5fb). This map reveals three different 
modes of current-induced dynamics. For small current 
amplitudes, in-plane precession (IPP) around the initial 
stable position is observed. The precessional angle rises 
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FIG. 5: (Color online) Dynamical phase diagram for 
Co(20)/Cu(10)/Py(4) /Cu(4) /Co(10)/IrMn(8) exchange- 
biased double spin valve as a function of the angle Q and 
normalized current density I/Io (with 7o = 10 8 Acm~ 2 ): (a) 
average value of the s z spin component; (b) dispersion of 
the s z spin component; (c) typical precessional orbits; (d) 
resistance oscillations associated with the IPP for SI = 0.47T 
and 7 = -1.37 (left part), and with the OPP for Q = O.Itt 
and 7 = — 1.37o (right part). The other parameters as in 
Fig.3. 



with increasing current amplitude. Above a certain criti- 
cal value of J, the precessions turn to out-of-plane preces- 
sions (OPPs). In a certain range of O, the OPPs collapse 
to a static state (SS), where the spin s remains in an 
out-of-plane position close to ±e x . 

As Sr, departs from the collinear orientation, the crit- 
ical current needed for destabilization of the initial state 
increases. This growth is mostly pronounced close to 

= 7r/2. To describe the critical current, we analyzed 
Eq. §E§ with respect to the stability of s in the upper po- 
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sition. Assuming that even in noncollincar configuration 
the stable position of s is close to 9 = 0, we have lin- 
earized Eq. |5]) around this point for arbitrary tt. Then, 
for the critical current needed for destabilization of the 
considered stable state we find 



i c,EBDSV 



afx M B d (ff ani + g "-+ g *» - h, 



a£(fi) - a^(Q)cosf2 



, (io) 



where a°(^) an( l °r(^) are calculated (for each config- 
uration) assuming 9 — > 0. Comparison of Eq. (|10|) with 
the results of numerical simulations is shown in Fig. [5] 
When considering the opposite stable point as the initial 
state, we need to take ol and or for 9 — ► it. Clearly, to 
destabilize the 9 — 7r state one needs current of opposite 
direction. 

Current-induced oscillations are usually observed ex- 
perimentally via the magnetoresistance effect. When 
electric current is constant, then magnetic oscillations 
cause the corresponding resistance oscillations, which in 
turn lead to voltage oscillations. The later are measured 
directly in experiments. In Fig. Hfd) we show oscilla- 
tions in the system resistance associated with the IPP 
(left) and OPP (right). As the amplitude of the oscilla- 
tions corresponding to the OPP mode is sufficiently large 
to be measured experimentally, the amplitude associated 
with the IPP mode is relatively small. This is the reason, 
why IPP mode is usually not seen in experiments. 



DISCUSSION AND CONCLUSIONS 



ing improvement has been found to be dependent on the 
spin-flip lengths in the magnetic and nonmagnetic lay- 
ers. According to our numerical simulations, an ultrafast 
subnanosccond current-induced switching processes can 
occur in DSVs with antiparallcl magnetic moments of the 
outermost magnetic films. 

In exchange-biased spin valves we have identified con- 
ditions which can lead to various types of spin dynam- 
ics. For fl < 7r/2, the negative current excites the central 
magnetic layer, while for f2 > tt/2, opposite current direc- 
tion is needed. We have evaluated parameters for which 
switching to a new stable state or to a precessional regime 
appears. This is especially interesting from the applica- 
tion point of view. However, the static state (SS) should 
be treated carefully. Such a state is often connected with 
the wavy-like angular dependence of ST T 13 i 16 i 17 . The 
SS becomes stable in the framework of macrospin model, 
when STT disappears in a certain noncollinear configura- 
tion. However, experimentally such a state has not been 
observed so far—. It has been shown in more realistic 
micromagnetic approach, that this static state may cor- 
respond to out-of-plane precessions^. The reason of this 
is the finite strength of exchange coupling in magnetic 
film, which docs not fully comply with the macrospin 
approximation in some cases. 

Finally, we have derived an approximate formula for 
the critical current in exchange biased DSVs, valid for 
non-collinear magnetic configurations. This formula, 
Ea.(|10p. is in good agreement with the numerical sim- 
ulations. 



We have calculated STT in metallic dual spin valves for 
arbitrary magnetic configuration of the system, but with 
magnetic moments of the outer magnetic films fixed in 
their planes either by large coercieve fields or by exchange 
anisotropy. In the case of symmetric DSV structures, we 
found a considerable enhancement of STT in the antipar- 
allel magnetic configuration. This torque enhancement 
leads to reduction of the critical current for switching as 
well as to reduction of the switching time. The switch- 
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